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General disagreement between the Geometrical Description of Dynamical In-stability 
-using non afRne parameterizations- and traditional Tangent Dynamics 
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In this paper, the general disagreement of the geometrical lyapunov exponent with lyapunov 
exponent from tangent dynamics is addressed. It is shown in a quite general way that the vector 
field of geodesic spread ^q is not equivalent to the tangent dynamics vector ^^ if the parameterization 
is not afSne and that results regarding dynamical stability obtained in the geometrical framework 
can differ qualitatively from those in the tangent dynamics. It is also proved in a general way that 
in the case of Jacobi metric -frequently used non afflne parameterization-, ^q satisfies differential 
equations which differ from the equations of the tangent dynamics in terms that produce parametric 
resonance, therefore, positive exponents for systems in stable regimes. 

PACS numbers: 



The definition of an appropriate indicator of dynamical 
(in)-stability for systems described by the hamiltonian 
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has been the subject of many works during the last cen- 
tury. We can cite, those primary methods based on per- 
turbation theory that appeared in the works of Poincare, 
usage of action variables, homoclinic intersections and 
qualitative analysis from the poincare surfaces sections 
PSS; or more quantitative and non-perturbative meth- 
ods like the Lyapunov exponent [l|, 0, S lj|- The latter, 
is a direct measure of the exponential divergence of close 
trajectories. It is computed as the limit 
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|q-(0,T2)-q(0,Ti) 

where r is a set of parameters which determines the initial 
conditions for each trajectory^. In this way A quantifies 
the exponent of divergence of two trajectories that have 
initial conditions ti and T2 close to each other. This in- 
dicator, although very intuitive, clearly fails for bounded 
systems, where trajectories can evolve very different but 
the distance ||q(t, r2) — q(i, ri)|| will never be greater than 
the size of the system. An alternative that was proposed 
to solve this problem was the tangent dynamics. This 
method, calculate a similar exponent but based on a dif- 
ferent measure. The distance ||q(i,T2) — q(t, ri)|| is re- 
placed by the norm of the solution Aq of the linearized 
variation of the equations of motion 5|]. 
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It is equivalent to follow the evolution of 
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a directional derivative which is taken at constant t^. 

Although the system might be bounded, the tangent 
dynamics accumulate the local gradients of divergence 
along the trajectory and it has proved to give a good 
measure of how sensitive a trajectory is with respect to 
variation of initial conditions. 

In the first half of the 20th Eisenhart, inspired by 
general relativity theory, proposed a geometrical descrip- 
tion of classical dynamics for systems which can be de- 
scribed by the least action principle Q. Geometriza- 
tion of dynamics represents the motion (actual trajec- 
tories) as a geodesic on a manifold with a suitable metric 
g = g^j^^dq^ ® dq" . For this purpose he introduced a met- 
ric that has been known after him as Eisenhart metric, 
whose arc length is given by 
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Here q^ = t and, from the integration of geodesies equa- 
tion, the extra coordinate is 
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Physical motions satisfy an affine parametrization ds^ — 
K^dt^ {k is a real arbitrary constant). 

It has been proposed that from a geometrical frame- 
work, stability analysis are also possible, by means of the 
geometrical lyapunov exponent q. 



l^Gis) 



Xg = lim -ln(- 

s^oo s ||€g(0)|| 



(7) 



'Corresponding author: |cuervo@inorg.chem.ethz.ch I 

^ In order to simplify the notation, here and in the following I omit 

the evident dependence of the exponent on the initial conditions 

and I write A instead of A(t) 



^ Directional in the sense that the parameter space has more that 
one dimension. Generally, if the system has N degrees of free- 
dom, 2N parameters are needed to determine the trajectory. 
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and it is called the vector field of geodesic spread. It is a 
directional derivative of geodesies at fixed arc length and 
evolves according to the Jacobi-Levi-Civita equation 
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where J| and ffi-^^ stand respectively for the covariant 
derivative along geodesies and the Riemann-Christoffel 
curvature tensor 8] . Eq. ([5]) , after opening the covariant 
derivative and the curvature tensor reads 



2efc 



ds2 



2rf,^ 



dq^ d^^ 



ij 



ds ds 



dq^ dq" 
' ds ds 



■e 



(10) 



It has been 



where Vf^ = ig''™ {dim,] + 9mj,i - 9ij,m)- 
shown that for the Eisenhart metric the spacial compo- 
nents of JLC equations are equivalent to the equations 
of tangent dynamics. 

Several works in the last decade have been dedicated 
to show that another metric can be used, in particular, 
the kinetic energy metric (also called Jacobi metric) @] 
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By applying the two mentioned metrics to some un- 
stable systems it has been argued the equivalence of 
the results @, [13], although in some cases a strange 
suppression of chaos with increasing number of non- 
separable degrees of freedom has been observed using 
Jacobi metric[9||. The fact that for many hamiltonian 
systems the curvature in the Jacobi framework is mostly 
positive has been used to support the idea that the main 
source of instability is the parametric resonance rather 
than the hiperbolicity of the potential surface^, M, \lM- 
There are some reasons already declared in previous 
works pointing at the inadequacy of the Jacobi metric. 
One reason is that it becomes singular at the bound- 
aries where the kinetic energy is zero. Thus, there are 
geodesies on the boundaries that do not correspond to 
physical motionsflj, [l3]. However, some authors con- 
sider that for many degrees of freedom, where the prob- 
ability of reaching the boundaries is practically zero, 
this metric gives good physical results, but it has not 
been proved. Some stimulating results regarding ge- 
ometrization in the thermodynamic limit have been ob- 
tained just within Eisenhart metric [ij, ll5|. A general 
proof of the equivalence of the geometrical approach and 
the tangen t dy namics, is still missing; while Cuervo and 
MovassaghjTyl have already shown the appearance of a 
positive lyapunov exponent in a trivially stable system 
when Jacobi metric is used. 

Here I will show briefly that the geometrical measure 
([5]) generally does not give results equivalent to those 



from the tangent dynamics when parameterization is not 
affine and I will explicitly show how the jacobi metric 
fails. 

Let g* be the coordinates used to describe a system, t 
the time and r the set of parameters that determine the 
initial conditions. The motions (trajectories) (/'(t, r) will 
depend of time and r. The arc length is a functional of 
the path; on extremal paths the arc length is a function 
of time and r. On geodesies, coordinates and arc length 
variations are 
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where {^\ is the variation of the arc length when one 
moves from one trajectory to another, at fixed time. 
(^)t is just a short form of (n- Vi-)t, a directional deriva- 
tive on the initial conditions taken at fixed time. From 
eqs. (fT2a|) and (J12bp 
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and using the definitions Eqs. @ and ([5]), one arrives to 
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Eq. p4)) relates in a compact way, the definitions of the 
divergence vector field obtained from tangent dynamics 
with that obtained from any geometrical approach. 

A consequence of Eg. p^ is that geometrical ^g and 
tangent dynamics ^t measures will differ unless \jf-\ 
vanish. This is the case of affine parameterizations, where 
the arc length (ds^ = K^dt^) depends just of the time 
elapsed. This is consistent with the fact that in Eisenhart 
metric k is an arbitrary positive constant independent of 
the trajectory and results regarding stability are always 
equivalent to those obtained from the tangent dynamics. 
In general cases, the directional derivative (ff)^ might 
have singularities or be ill-defined. 

It is evident that if (ff ) has zeros along the trajectory 
the geometrical measure will diverge. This is the case of 
Jacobi metric ((§f) = 2r) when the kinetic energy can 
be zero (very common in low dimensional systems). How- 
ever, systems with more degrees of freedom, where zeros 
in the kinetic energy are not likely to happen, require 
a further analysis. For it, one can look at the differen- 
tial equations satisfied by and ^^. For the Tangent 
dynamics the equations of motion are 
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while for the Jacobi measure, after substituting F^ = t^ + y {a''™'aijVm — S'-Vi — SfVj) in pU|) one obtains 
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The equation for ^j contains the same terms that the one 
for ^T plus other terms which have been written in the 
second hne. The first conclusion that can be drown from 
it is the already known appearance of divergent terms 
when kinetic energy is zero. In addition these terms, as 
will be shown immediately, introduce non physical insta- 
bilities also for trajectories where kinetic energy never 
vanishes. 

Many of the systems of interest are bounded. The 
phase space has a finite volume. For energies close to a 
minimum of the potential, the motion will be restricted 
to small oscillations around an equilibrium point. In this 
case, harmonic approximation for the potential energy is 
valid and systems are dynamically stable with periodic 
trajectories. Since the stability analysis should not de- 
pend of the variables used, let expand for simplicity, the 
potential around a minimum in cartecian coordinates. 
So, aij are constant, i^ = 0. 

Tangent dynamics equations become 

C^ = -a^'^Vkier (17) 

^ Expansion of V up to n-order (n > 2) in q gives V ^i of order 
n- 2 



For up to second order expansion in the potential V{q) 
these equations give oscillatory solutions for £^t leading 
to vanishing exponent. When the system is more excited 
above the minimum energy, terms of order higher than 
two in coordinate dependence of the potential must be 
considered. Here motion will be quasi-periodic and this 
regime corresponds to the onset of chaos. Vu (the fre- 
quencies of the tangent dynamics) will be q-dependent 
(time dependent, periodic)^. Since Vki oscillate at har- 
monics of the frequencies of the system they produce 
parametric resonance in the tangent dynamics and there- 
fore positive exponents. For very high energies, systems 
start filling the phase space, trajectories are no longer 
periodic but unstable due to hyperbolic points between 
the minimums of the potential. Thus, parametric reso- 
nance in ^ dynamics is relevant at the on-set of chaos in 
hamiltonian systems. 

The picture is very different in the Jacobi framework. 
The extra terms in the equations of motion for ^j, (see 
bellow in cartecian coordinates eq.((T8])) 
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contain q-dependent terms also when the potential is 
quadratic in q, which therefore oscillates with the fre- 
quencies of the system and its harmonics. For even po- 
tentials (it is always the case in the harmonic approxima- 
tion), the kinetic energy oscillates with frequencies that 
are twice the frequencies of the system. So, the extra 
terms in eg. lfTS)) set the perfect scenario for parametric 
resonance, and it has nothing to do with non-linear be- 
havior of the physical system. For systems at the onset of 
chaos, real parametric resonance appears from the first 
term of the rhs in Eq. (|18p and for higher energies the 



effect of the hyperbolic points will come into the game. 
Thus, it is not possible to discriminate the false exponen- 
tial divergence from the physical one, when the exponent 
is computed from the evolution of S,j- Stable regimens 
can be shown as chaotic. 

Another point to stress is that Jacobi and Eisenhart 
metric are not related by a canonical transformation, 
they are defined in different spaces. The statement 
"physics is independent of the frame of reference" , is then 
out of context since these coordinates must be defined on 
a same space and related by a canonical transformation. 



In addition, the independent variable (arc length) cho- 
sen in the Jacobi metric is not invariant with respect to 
canonical transformations. 

In conclusion, whether the physical system is chaotic or 
not, and the nature of this instabilities (hyperbolic points 
or parametric resonance) are questions that can not be 
answered by Jacobi metric. It always gives positive expo- 
nents for hamiltonian bounded systems. Other non-affine 



parameterizations might have similar drawbacks. In gen- 
eral, the vector field of geodesic spread and the tangent 
dynamic vector field are not equivalent when the param- 
eterization with the arc length is not affine. 
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